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Abstract 

The numerically dominant QED contributions to the anomalous magnetic mo¬ 
ment of the muon stem from Feynman diagrams with internal electron loops. We 
consider such corrections and present a calculation of the four-loop light-by-light- 
type corrections where the external photon couples to a closed electron or muon 
loop. We perform an asymptotic expansion in the ratio of electron and muon mass 
and reduce the resulting integrals to master integrals which we evaluate using an¬ 
alytical and numerical methods. We confirm the results present in the literature 
which are based on different computational methods. 

PACS numbers: 12.20.-m 12.38.Bx 14.60.Ef 


1 Introduction 


The anomalous magnetic moment of the muon provides an important test of the Stan¬ 
dard Model of particle physics. It has been measured to an impressive accuracy at the 
Brookhaven National Laboratory mca and two new experiments at Fermilab ||3] and 
J-PARC [l] are planned to further improve the measured value. 

On the theory side much effort has been made to provide a precise prediction for the muon 
magnetic moment, see e.g. Refs. EH3 for detailed reviews. The comparison of the precise 
measurements and calculations shows a deviation of about three standard deviations, 
which already persists for several years. This fact makes the anomalous magnetic moment 
of the muon, a M , an interesting quantity for further investigations. 

Several ingredients are needed to obtain the theory prediction for a M . The numerically 
most important one origins from QED radiative corrections which are analytically known 
up to three loops IMS] and numerically up to five-loop order du Also the electroweak 
correction, which are known at the two-loop level, are under control P2H3S1. The dominant 
contribution to the uncertainty comes from the hadronic contribution which can be split 
into a vacuum polarization and light-by-light contribution. The vacuum polarization 
contribution is obtained with the help of a dispersion integral over the experimentally 
measured cross section e + e _ —> hadrons where the dominant contribution comes from 
low energies. The corresponding analysis has been performed at leading order mm, 
next-to-leading order PHE0H22] and next-to-next-to-leading order [23j. The least known 
contribution origins from the hadronic light-by-light part which has been considered by 
several groups at leading order [24H26] . The corresponding next-to-leading order effects 
have been estimated to be small [27 ], 

In this paper we focus on the QED contribution to the muon anomalous magnetic moment 
which can be cast in the form 



n =1 


(i) 


where a is the fine structure constant. The first three coefficients on the right-hand- 
side, which correspond to the one-, two- and three-loop corrections, are known analyti¬ 
cally |8l ll0ll28l43T] . For the four- and five-loop contributions only numerical results are 
available [ 1111351136 ] . Note, that even for the four-loop coefficient a ® there is no systematic 
cross check by an independent calculation; only a few special cases have been computed 
analytically (see, e.g., Refs. [3ZH2QJ). An independent calculation of a® is important 
since the four-loop contribution in Eq. (pQ) amounts tc@ 
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The ellipses stand for further digits and small contributions which are not shown. 
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Figure 1: Sample light-by-light-type Feynman diagrams contributing to a M . The external 
solid line represents the muon and at least one of the internal solid loops denotes electrons. 
In the case of IV(a) the second fermion loop can either be an electron or muon loop. Wavy 
lines represent photons. 

which is comparable to the deviation between the experimentally measured and theoret¬ 
ically predicted result for given by HU 

a M (exp) — a^(SM) « 249(87) x 10 -11 . (3) 

In Eq. ([2]) we have separated the contributions containing at least one closed electron 
loop (second term on left-hand-side) from the pure photonic part. The former are nu¬ 
merically dominant^ which provides the motivation to concentrate in a first step on these 
contributions. Note that about 95% of the electron loop contribution originates from the 
so-called light-by-light-type Feynman diagrams where the external photon couples to a 
closed fermion loop. Such contributions arise for the first time at three loops; see Fig. |T] 
for sample diagrams. In this paper we perform an independent calculation of the four-loop 
corrections. 

It is convenient to decompose the four-loop term a/ ; into a purely photonic piece and 
contributions involving electron and/or r loops. Following Ref. [TT] we write 

a (s } = 4 ( 8 ) + 4 ( 8 ) (m^/m e ) + Af ] {rnJm T ) + A® (mjm e , m M /m r ). (4) 

^ 2 8 \ rn n/m T ) has been computed in Ref. [41] using an asymptotic expansion for m 2 <C 
rri 2 . Analytic results have been obtained for several expansion terms which show a rapid 
convergence. A^^m^/m-r) and Af\rn fl /m e ,m li /m T ) are suppressed by m 2 ^/m 2 and thus 
they are numerically small. 

A^\m,j,/m e ) can be split into light-by-light-type contributions (cf. Fig. [T]) and contribu¬ 
tions where the external photon couples to the external muon line. The leading term of 
the latter can be obtained from calculations where in a first step the electron mass is set 
to zero and the fine structure constant is renormalized in the MS scheme. Afterwards 
a is transformed to the on-shell scheme which introduces log(m e /m M ) terms in the final 

2 This is also true at two and three loops, see. e.g., Ref. [6] for explicit results. 
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result. Using this approach, the non-light-by-light contributions with two closed electron 
loops have been computed analytically in Ref. [39]. 

/o\ 

In this paper we compute the four-loop light-by-light contributions to HI, (m^/rrie) which 
are exemplified by three Feynman diagrams in Fig. [I] In case the external photon couples 
to a closed electron loop it is not possible to set m e = 0 since this generates infrared 
singularities. To circumvent this problem we perform an asymptotic expansion for m e <C 
which is described in some detail in Section [2] Results for contribution IV(a) with 
two closed electron loops have already been considered 40 years ago in Refs. [32133] • In 
Section |3] we will discuss in detail our results and compare to their findings and also to 
the ones in Ref. m Section [4] contains our conclusions. 


2 Technical details 

/o\ 

The Feynman integrals which contribute to the light-by-light part of H), {m^/m e ) contain 
two widely separated scales, which provides a small expansion parameter 

Tfl 

x = — « 1/206.7682843. (5) 

m^ 

Thus, it can be expected that already a few expansion terms provide a good approximation 
to the exact result. We compute four terms and show that the one of order x 3 leads to 
negligible contributions. The linear and quadratic term, however, can still lead to sizable 
contributions since the coefficients of x n contain log(x) ~ —5.3 terms which at four-loop 
order are raised up to forth power. 

We have implemented the asymptotic expansion using two different programs. In the first 
approach we use the Mathematica package asy [3US5] which is based on expansion by 
regions [46]i47] formulated at the level of the alpha representation [48]. It provides the 
possibility to obtain the asymptotic behaviour of a Feynman diagram in a given limit. In 
fact, the output of asy are scaling rules for the alpha parameters. To exploit this infor¬ 
mation one has to find a distribution of the external momentum and the loop momenta 
obeying certain scaling rules which we obtain by trying out all possible combinations. 

The second approach is based on an in-house program which generates all possible com¬ 
binations of loop momenta and external momentum and assigns for each combination 
all possible scalings of the loop momenta (i.e. each loop momentum can either be soft 
or hard). In this way one obtains by construction all contributing regions. However, a 
double counting is introduced since the routing of the loop momenta is not unique. The 
double counting is eliminated with the help of the unique alpha representation which is 
generated for each momentum distribution. 

We have applied both methods to each Feynman diagram and have obtained identical 
final results. Note, however, that asy requires significantly more CPU time than our 
in-house program which is tailored for the problem at hand. 
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Figure 2: Sample Feynman diagrams for the three appearing integral classes introduced 
in Eqs. (JZJ) (left), (J8]) (middle) and (EJ) (right). Solid and dashed lines denote massive and 
massless propagators. Double lines stand for denominators of the form 1/(2 k ■ q ). 


The output of the asymptotic expansion is manipulated with FORM [J9]l5l3] and TFORM j5T] 
(see also Ref. [52]) which are used to perform traces and to deal with tensor structures. 
Afterwards only scalar integrals are left which are reduced to master integrals with the 
help of FIRE [53] and crusher [54]. For the classification of the master integrals we 
introduce the following set of propagators 


Pi=<To, 


k 2 ' k 2 -M 2 ’ k 2 - 2k- q 2k- q 


M G {m e , m^} , 


( 6 ) 


where k is a linear combination of loop momenta and g, the external momentum, with 
q- = m 2 . Using these propagators we can build the required three integral classes 


Vacuum integrals 


U(ii,..., in') 



Di G {Pi, P2} , 


(7) 


On-shell integrals 

■ ■., in) 





Di e {Pi, P 2 , P3}, (8) 


Linear integrals 

P(U) • • • An) 





Di G {Pi, P 2 , P4} , 


(9) 


which are exemplified in Fig. [2] 

To get an impression how the individual types of integrals arise after asymptotic expan¬ 
sion we discuss in some detail the three-loop case (cf. left diagram in Fig. [[]). The valid 
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regions are obtained by considering appropriate routings of the loop and external mo¬ 
menta, allowing each loop momentum to be either soft (k ~ m e ) or hard ( k ~ m^). In 
total we obtain eight possible regions. In case all loop momenta are hard the electron 
propagators are expanded in m e and one ends up with three-loop on-shell integrals. On 
the other hand, in case all loop momenta are soft the muon propagators are expanded for 
k 2 <C 2 k ■ q and one has 

1 1 1 / k 2 

m ‘fi ~ ~~ <?) 2 2 k ■ q — k 2 2 k ■ q \ 2 k ■ q + 

which leads to linear integrals where m e sets the mass scale. If two loop momenta are hard 
and one is soft one obtains two-loop on-shell (with mass scale m and one-loop vacuum 
integrals (with mass scale m e ). The remaining three regions, i.e. two soft and one hard 
loop momentum, leads to one-loop on-shell integrals and either vacuum or linear integrals 
where the massive scale is given by the electron mass. It is interesting to note that on- 
shell and vacuum integrals only occur for the even powers of x whereas linear integrals 
are present both for even and odd powers. Note, however, that the corresponding master 
integrals differ: the x° and x 2 terms involve master integrals with even number of linear 
propagators whereas for the x l and x 3 terms their number is odd. 

The pattern observed at three loops repeats itself at four loops. For some diagrams one 
obtains more than 20 regions leading to single-scale integrals of the type 

• four-loop on-shell, 

• four-loop linear, 

• products of three-loop on-shell and one-loop vacuum, 

• products of two-loop on-shell and two-loop vacuum or linear, 

• products of one-loop on-shell and three-loop vacuum or linear, 

• and products two one-loop and one two-loop integrals involving vacuum, on-shell 
and/or linear integrals. 

The mass scale of the on-shell integrals is given by and the one of the vacuum and 
linear integrals by m e , with the following exception: one-loop vacuum integrals with mass 
scale m^ occur in IV(a) in case a closed muon loop is present. 

Massive vacuum integrals, which are only needed up to three loops, are well documented 
in the literature [55}l56j. All one-, two- and three-loop on-shell integrals entering our 
calculation are known analytically 1571 • There are about 70 four-loop on-shell integral 
which are needed for IV(a), IV(b) and IV(c), about 40 are known analytically or to high 
numerical precision. All of them are taken from the calculation of the four-loop MS- 
on-shell quark mass relation performed in Ref. |58j . We furthermore require about 70 
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four-loop linear master integrals where about 20 have been computed analytically or to 
high numerical precision using Mellin-Barnes techniques [59] . The remaining ones have 
been computed using the package FIESTA [60] • We propagate the uncertainty from each 
e coefficient of each master integral to the final result. To obtain the final error estimate 
we add the uncertainties in quadrature. 

The renormalization of the light-by-light contribution only involves one-loop on-shell coun¬ 
terterms for the fine structure constant, the muon and electron masses and the wave 
function which are well established in the literature (see, e.g., the textbook my 


3 Discussion of results 


We start with looking at the three-loop light-by-light contribution to A^irn^/me) which 
is known analytically [8]. Using the method described in the previous section we obtain 
in numerical form (x = m e /m 


A 


( 6 ) 

2,lbl 


( m^/rrie ) 


-14.5525 - 6.57974 
+ x [18.0169- 13.15954] 

+ x 2 [-12.4582 + 1.59444 - 5.5700^ + 0.6667^] 

+ x 3 [-12.0628- 10.96624] 

+ x 4 [14.4529 - 19.04314 + 2 . 2778£ 2 - 0.7778^] 

-14.5525 + 35.0805 
+ [0.08714 + 0.3393] 

+ [-0.0002914 - 0.0001988 - 0.003703 - 0.002363] 

+ [-0.0000014 + 0.0000066] 

+ [0.0000000079 + 0.000000056 + 0.000000035 + 0.000000064] 
20.5280 + 0.4265 - 0.006557 + 0.0000052 + 0.0000002 
20.9479, (11) 


where terms of order x 5 are neglected and 4 — log(:r) ~ —5.3. For completeness we 
provide analytic results in the Appendix. One observes that odd powers in x contain at 
most linear terms in 4 whereas even powers contain terms up to £ 3 . After the second 
equality sign the numerical value for x is inserted, however, the contribution from the 
various powers of x and 4 are kept separately. At order x° and x l the logarithmic 
contribution dominates over the constant. At order x 2 the constant and linear logarithmic 
term is of the same order of magnitude and about a factor ten smaller than the quadratic 
and cubic contribution. After the third equality sign all contributions to x n are added. 
One observes that the overall contribution of the x 2 term is already quite small and 
amounts to only 0.03% of the leading term. The linear term provides a 2% contribution 
and is still important. Let us mention that the cubic and quartic terms are below 0.00003% 
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and are thus negligible. For completeness we present the final result for A^ bl (m /t /m e ) 
after the last equality sign. 

We now turn to the four-loop results. For convenience we split IV(a) into three contribu¬ 
tions: IV(aO) contains two closed electron loops, in IV(al) only the fermion loop with the 
coupling to the external photon contains electrons, and in IV(a2) electrons are running 
only in the two-point polarization function and muons in the other fermion loop. Let us 
mention that the coefficients of the logarithmic contributions in the case of IV(aO) are 
known analytically since the contributing four-loop on-shell integrals are available in the 
literature [32] and only the four-loop linear integrals have to be evaluated numerically. 
As a consequence, it is possible to reconstruct the pole terms of the IV(aO) contribution, 
which are in one-to-one correspondence to the coefficients of the logarithms, analytically. 
Similar arguments can be used to obtain the logarithmic contributions of IV(al) and 
IV(a2), for the x 1 and x 3 terms of IV(b) and for the a: 0 , x l and x 3 terms of IV(c). They 
are given in the Appendix. In the main part of the paper we restrict ourselves to numerical 
results. 

Using the same scheme for presenting the results as in Eq. (HID we obtain for the five 
four-loop light-by-light-type contributions 


^(8),IV(aO) 


^(8),IV(al) 


7.5018 ± 0.0026 + 14.88084 + 6.5797^ 

+ x [6.29 ± 0.46 - 14.62164 + 8.7729^] 

+ x 2 [-16.81 ± 0.43 + 30.01724 - 6.5069^ + 7.6489£ 3 - 0.8889^] 

+ x 3 [-48.31 ± 0.24 - 4.87394 + 13.1595^] 

7.5018 ± 0.0026 - 79.3384 + 187.0352 
+ [0.0304 ± 0.0022 + 0.3770 + 1.2061] 

+ [-0.000393 ± 0.000010 - 0.003743 - 0.004326 - 0.02711 - 0.01680] 

+ [-0.0000055 + 0.0000029 + 0.0000423] 

[115.1986 ± 0.0026] + [1.6135 ± 0.0022] + [-0.052378 ± 0.000010] 

+ [0.000040] 

116.7598 ±0.0034, (12) 

2.734 ±0.028 
±x [-9.5571] 

± x 2 [-7.494 ± 0.033 - 14.20104 ± 0.3559£ 2 ] 

±x 3 [6.1668 -9.06544] 

2.734 ±0.028 
± [-0.04622] 

± [-0.000175 ± 0.000001 ± 0.001771 ± 0.0002366] 

± [0.000001 ± 0.000005] 

[2.734 ± 0.028] ± [-0.04622] ± [0.001832 ± 0.000001] ± [0.000006] 


(13) 


^(8),IV(a2) 


^(8),IV(b) 


^(8),IV(c) 


2.690 ±0.028, 

0.370 ± 0.033 - 0.74204 
± x [ 0 ] 

+ x 2 [0.8628 ± 0.002618] 

± x 3 [-2.6844] 

± [0.370 ± 0.033 ± 3.9561] 

+ [ 0 ] 

± [0.00002018 ± 0.00000006] 

± [-0.0000003] 

± [4.326 ± 0.03288] ± [0] ± [0.00002018 ± 0.00000006] + [-0.0000003] 

4.326 ±0.033, (14) 

27.395 ± 0.014 ± (4.93482 ± 0.00003)4 
± x [-0.81 ± 1.22 ± 59.02354] 

± x 2 [142.5 ± 7.6 ± 40.65464 ± 20.55824 - 9.61674 ± 0.83334] 

± x 3 [62.11 ± 2.89 ± 132.74214 - 40.94064] 

27.395 ± 0.014 ± (-26.3105 ± 0.0002) 

± [-0.0039 ± 0.0059 - 1.5219] 

± [0.003334 ± 0.0001769 - 0.005070 ± 0.01367 ± 0.03409 ± 0.01575] 

± [0.000007 ± 0. - 0.000080 - 0.000132] 

[1.084 ± 0.014] ± [-1.5259 ± 0.0059] ± [0.06177 ± 0.00018] ± [-0.0002047] 
-0.380 ±0.016, (15) 

-14.900 ± 0.059 - 3.28994 
± x [65.4209] 

± x 2 [33.61 ± 9.96 ± (56.76 ± 0.78)4 ± (-24.66 ± 0.27)4 
±(-1.935 ±0.091)4] 

± x 3 [25.11 ± 0.60 ± 107.43254 - 19.86104] 

-14.9 ±0.05945 ±17.54 
± [0.3164] 

± [0.00079 ± 0.00023 ± (-0.007078 ± 0.000097) 

±(-0.01640 ± 0.00018) ± (0.00686 ± 0.00032)] 

± [0.00000284 ± 0.00000007 - 0.00006479 - 0.00006387] 

[2.641 ± 0.059] ± [0.3164] ± [-0.01583 ± 0.00045] 

± [-0.00012582 ± 0.00000007] 

2.941 ±0.059, (16) 


where the uncertainties origin from the numerical integration using FIESTA. It is common 
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A^im^/me) 

this work 

HUES] 

m 

03 

IV(aO) 

116.76 ±0.02 

116.759183 ±0.000292 

111.1 ±8.1 

117.4 ±0.5 

IV(al) 

2.69 ±0.14 

2.697443 ± 0.000142 



IV(a2) 

4.33 ±0.17 

4.328885 ±0.000293 



IV(a) 

123.78 ±0.22 

123.78551 ±0.00044 



IV (b) 

-0.38 ±0.08 

-0.4170 ±0.0037 



IV (c) 

2.94 ±0.30 

2.9072 ± 0.0044 




Table 1: Summary of the final results for the individual four-loop light-by-light-type 
contributions and their comparison with results from the literature. Note that the uncer¬ 
tainties given in the second column are obtained from Eqs. da to (11611 after multiplication 
by five. 


to all five cases that the x 3 term only provides a negligible contribution of at most 0.05% 
[in the case of IV(b)] which is much smaller than the uncertainty estimate of the leading 
x° term. In most cases the x 2 terms lead to contributions comparable to the numerical 
uncertainty of the leading term. Together with the x 3 terms they confirm the good 
convergence property of the asymptotic expansion. It is nevertheless interesting to note 
that in general the cubic and quartic logarithms of the x 2 term lead to the largest numerical 
contributions, however, also the quadratic log(x) terms are not negligible. 

In the linear mass correction terms we observe logarithmic contributions up to second 
order. They lead to significantly larger numerical contributions than the linear logarithms 
which are in turn much larger than the constant. 

In the leading term only IV(aO) has a quadratic logarithm, one from each electron loop. 
The log 2 (±) term provides the numerical dominant contribution, however, more than 40% 
are canceled by the linear logarithm; the constant term is an order of magnitude smaller. 
Linear logarithmic terms are also present for IV(a2), IV(b) and IV(c). In the latter two 
cases strong cancellations between the log(x) and the constant are observed, which leads 
to an interesting effect for IV(b): the leading term is smaller than the mg/m^-suppressed 
term. As a consequence the final numerical result for IV(b) is quite small and has a big 
relative uncertainty. 

In Table |T] we summarize our Endings for IV(aO), IV(al), IV(a2), IV(b) and IV(c) and 
compare with the literature. To be on the conservative side we multiply the Monte-Carlo 
uncertainty from Eqs. (1T21) to (fTfill by a factor of five. For all classes under consideration 
we find good agreement with previous results. For the dominant contribution IV(a) we 
find excellent agreement. 

Although our numerical precision can not compete with the one of Ref. mi, let us note 
that our computation procedure is completely different from the numerical method used 
in HU. As our calculation is a second evaluation of the complete light-by-light contribution 
to A^\m li /m e ), our result is an important check of the existing value. Furthermore, since 
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the asymptotic expansion and the reduction to master integrals were done analytically, it 
is possible to systematically improve the precision by evaluating more and more master 
integrals analytically. 


4 Conclusions 

We have computed the four-loop QED corrections to the muon magnetic moment induced 
by the light-by-light-type Feynman diagrams involving a closed electron loop. This is 
a finite and gauge invariant subset, which provides about 95% of the total four-loop 
contribution. Our results for the individual sub-classes IV(a), IV(b) and IV(c) agree with 
the literature HU- 

We want to stress that our approach is completely different from the one of Ref. HU 
and thus should be considered as an independent cross check. We perform an analytic 
reduction of all occurring integrals to a small set of master integrals which are then com¬ 
puted using analytical or numerical methods. In particular all counterterm contributions 
are available analytically. On the other hand, in Ref. HU infrared and ultraviolet finite 

(o\ 

multi-dimensional integrands for the individual sub-classes of ajj ,' are constructed which 
are then integrated numerically. 

Our final results can be found in Table [U Although our numerical uncertainty, which 
amounts to approximately 0.4 x (a/7 r) 4 1.2 x 10 —11 , is larger than the one of Ref. HU, it 
is still much smaller than the difference between the theory prediction and experimental 
result of [cf. Eq. (J3])]. We also want to note that our result can be systematically 
improved by evaluating more and more master integrals analytically. 
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Appendix: Selected analytic three- and four-loop re¬ 
sults for an 


In this appendix we provide analytic results for the three-loop light-by-light contribution 

A' 


fibiK/ m e) and f° r some °f the four-loop contributions. 


The analytic version of Eq. (1TT1) reads 

107T 2 597T 4 




X 


/ 424n 2 


+ X 


283 


270 

196 ln(2)7r 2 
3 

25vr 2 


3Cs 


27r 2 


47T 2 


61 


12 

32t r 2 


+ 


18 

167T 4 


6l7T 4 47T 2 (3 

w + 3C3 + — 


+ X 


X 


,3 

9 

+ 135 

3 

f 117r 2 

107T 2 


l 9 

9 


f13283 

191t r 2 

4 

—— + — 


+ 4Cs) 4+(y-^ 


20 

"3 


2592 


+ 1 + 


517 13tT 


108 


9 


fi 2 + - 

3 


41 


, j_ F - - 

x + 18 x 9 


(17) 


where is Riemann’s zeta function evaluated at n. 

At four-loop order one can obtain analytic results for the coefficients of (n > 1) for 
IV(aO), IV(al) and IV(a2), for the x 1 and x 3 terms of IV(b) and for the x°, x 1 and x 3 
terms of IV(c). These contributions are given by 


A 


(8),IV(aO) | 


A 


(8),IV(al) i 


4 87T 2 
3 ^ 3~~ 


+ x 


X 


+ X 


87T 4 

" 45 " 


2 tt 2 

+ 6 C 3 ) 4 H- 7T 


-407T 2 


27 

47 


+ 


8 tt 2 


9 


2 987T 4 87r 2 £ 3 

3 +6 ^-i5B- 10& -^r 


5oc 5 


220 32t r 2 

+ 


9 9 

—407r 2 


16tt 4 14Cs 
" l35~ 3~~ 


+ 


82 

~9~ 


47T 2 

~27 


= x 


81 

137 12629tt 2 

9720 


47T 2 

4+ 


t 


21n(2)7r 2 31 tt 4 125C 3 \ 

3 "405" + 27 ) 


P 3 _ _ t 

-x g C X 


(18) 
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(19) 


X 


71 - 


20 


C-x i 


124vr 2 


9 27 J x 135 

^( 8 ),IV(, 2 )| = I _ 10 _ 32 _ 4hA2) _ 931^ 2 

9 3 27 K J 

41n 2 (2)7T 2 4l7r 4 8(3 , 57T 2 C3 5(5 


A 


(8),IV(b) I 


si A, = ® 3 tt - 


3 270 

317T 2 ' 


3 + 9 3 


9 


A 


(8),IV(b), 




137tt 2 1153t r 3 \ „ 112 tt 2 - 

+ _ I 4 - 


4050 


270 


27 


A 


(8),IV(c) I 


71 


?°/x g 4 j 


A 


(8),IV(c) | 


\X 


= 0 , 


^(8),IV(c)| 




2939vr 2 

270 


163tt 2 

81 


with <24 = Li 4 (l/ 2 ). 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 
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